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Abstract 



The purpose of this paper is to analyze the singularities of a well known bench- 
^ . mark problem "Andrews' squeezing mechanism". We show that for physically rel- 

I evant parameter values this system admits singularities. The method is based on 

■ Grobner bases computations and ideal decomposition. It is algorithmic and can thus 
be applied to study constraint singularities which arise in more general situations. 
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1 Introduction 

X. 

■ The "Andrews' squeezing system" was first described by Giles in |Gil78] and further studied 
in |Man81j . It is a planar multibody system whose topology consists of closed kinematic 
loops (see Figure [1]). The Andrews' system was promoted in ^SchQOj as a benchmark 
problem to compare different multibody solvers. Nowadays it is a well-known benchmark 
problem |HW91t IMI03j for numerical integration of differential-algebraic equations as well. 
The equations are of the Lagrangian form (or descriptor form, see also |Arp01| ) 

'f{t,y,y',y",X)=0 

^9{y) = ^ ^ 

where the function / describes the dynamical equations and g gives the (holonomic) 
constraints. Here y G M" are the (generalized) position coordinates, y' and y" are the first 
and second derivatives, respectively, and A is the Lagrange multiplier. 
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It is well known that singularities of any kind hinder solving equations numerically 
|RS88t IHW911 IBA941 IEH95] . Intuitively, a singularity is where the (generic) number of 
degrees of freedom of the system changes. Mathematically these are the points where 
the rank of the Jacobian of g drops. Hence in this paper we will not consider the actual 
dynamical equations and analyse only the constraints given by g. 

Most differential equation solvers include a possibility to monitor singularities, and 
usually when proximity of a singularity is detected, the computation is best to be inter- 
rupted. But this kind of monitoring is local only, that is, it does not tell us a priori where 
the singularities lie but only alert us when it is too late to fix things, so to speak. Also, 
the monitoring is often a non-negligible part of computational cost. Therefore, it would 
be highly useful to know a priori where the singularities are, or to make sure that there 
are no singularities, or perhaps even remove them (for the latter approach, see |Arp01| ). 
Locating singularities has been studied also in [McCOOj . If we cannot avoid or remove the 
singularities, at least knowing where they are encountered is helpful (indeed, necessary) 
when planning the computation without interruptions. One can then tune the chosen 
integration algorithm such that the disturbing effect of the singularities is diminished, 
for example by compensating the singularity of the Kepler problem by a local change 
of variables as in |LR05] within the computation. Further techniques on compensating 
singularities in multibody systems are gathered and concisely compared in [BA94j and 



The paper is organized as follows: in the next Section we present the situation in 
detail and formulate the constraint equations in polynomial form. Section [3] gathers the 
necessary algebraic tools. Section H] contains the actual analysis where we show that the 
mechanism indeed has singularities for certain parameter values. In Section [5] there are 
some numerical examples of singular configurations, and in Section [6] we summarize and 
discuss the results, and address possible future work. 

2 Andrews' squeezing mechanism 

The squeezing mechanism is given by the following equations. 



ai sin(yi) - 02 sin(?/i + 2/2) + as cos(?/3) - 62 
ai cos(?/i) - a2 cos(?/i + ^2) - ^4 sin(?/4 + y^) - cos(?/5) - Wi 
ai sin(?/i) - 02 sin(?/i + 2/2) + ^4 cos{y^ + y^) - ag sin(?/5) - W2 
ai cos(?/i) - 02 cos(?/i + 112) - ae cos(?/6 + Vr) - a? sm^y^) - wi 
^ai sin(?/i) - 02 sin(?/i + 1/2) - sin(?/6 + Vr) + 07 cos(?/7) - W2 



Compared to the original articles mentioned above, we have chosen the following notation 
for the parameters and angles: 



pT95] . 




ai = rr a2 = d = ss 04 = e 0,5 = zt = zf a-j = u 
hi = xb 62 = yb wi = xa W2 = ya 

yi= P y2 = Q 1/3 = 7 1/4 = $ 2/5 = 5 2/6 = ^ 2/7 = e 



so the positions in Cartesian coordinates of the fixed nodes A and B are given by 6 = 
(61,62) and w = {wi,W2), and the lengths of the rods by a = (ai, . . . ,07), see Figures [U 
and El 



Fixing the parameters a, b, and w, we have a map g : M'' — > M^. Hence the set of 
possible configurations, which is the zeroset Mg = g~^{0), is in general a curve (or possibly 
empty). Our task is to analyse the singularities of Mg, so let us state more precisely what 
is meant by a singularity. As mentioned before, in a singularity the number of degrees 
of freedom changes. It is well known |RS88l IBA941 IMcCOO] that this corresponds to the 
situation where the rank of Jacobian drops. 

Definition 2.1. Let / : ^ be any smooth map where k < n and let df be its 
Jacobian matrix. Let M = /~^(0) C be the zeroset of /. A point g G M is a singular 
point of M, if df does not have maximal rank at q. 

What in fact geometrically "happens" at a singular point may be quite complicated 
to determine. Typically the tangent space to M does not change continuously in the 
neighbourhood of a singular point, or possibly M intersects itself there. However, in all 
cases numerical problems occur, so it is important to try to find all singular points. 

Note that the constraint equations ([2]) (and hence the elements of its Jacobian ma- 
trix) are not polynomials, yet our algebraic approach works only in a polynomial setting. 
However, this problem is circumvented by reformulating g{y) as polynomials in the sines 
and cosines of yi by using the trigonometric identities 

cos(x)^ + sin(x)^ = 1 

sin(a; ±y) = sin(a;) cos{y) ± cos(x) sm{y) 

cos{x ±y) = cos(a;) cos{y) =F sin(a;) sin(?/) 

Setting Cj = cos(?/j), Sj = sin(?/j) we get the equations 

fliCi — a2{ciC2 — S1S2) — cisSa — h = 
aisi - a2{siC2 + C1S2) + a^c-i - 62 = 

fllCi - 02(0102 - S1S2) - a4(s4C5 + C4S5) - OsCs - Wi = 
p{0, S) = { OiSi- 02 (siC2 + C1S2) + 04 (C4C5 - S4S5) - 0585 -W2 = (3) 
OiOi — 02 {01O2 — S1S2) — Oq {0Q0^ — SqSt) — 0787 —Wi=0 
OlSi — O2 (siC2 + C1S2) — Oq [sqOy + OqSj) + OjOy — W2 = 

c2 + - 1 = 0, i = 1, . . . , 7. 



We have 13 polynomial equations {pi = 0), 11 parameters (ai, ... ,07, bi,b2,wi,W2) and 
14 variables (ci, si, . . . , C7, S7). Note that each pi is of degree two in q, Sj. The equations 
Pi = 0,...,Pq = correspond directly to the 6 original equations g{y) = with the 
simple substitutions above (for example cos(?/i + 7/2) = C1C2 — S1S2) and the equations 
P7 = 0, . . . ,pi3 = are the extra identities due to "forgetting" the angle variables yi. 

Note that this reformulation of the constraints as algebraic equations is not just a 
trick which happens to work in this special case; indeed most constraints appearing in the 
simulation of multibody systems are of this type. 

Now the above equations define p as a map p : M^^ — R^'^. Hence we expect that the 
zeroset V = p~^{0) C R^^ is a curve (or possibly empty). Singularities are then the points 
of this curve where the rank of dp is not maximal. To find these points we need now to 
introduce some tools from commutative algebra. 




Figure 2: The lengths aj and nodes of the Andrews' system. 



3 Background 



In this section we present briefly the necessary definitions from commutative algebra and 
algebraic geometry. More details can be found in |CL092] . |GP02j . [Nor 76] . and |Eis96]. 
These are roughly in the order of increasing difficulty, [CL092j being the most accessible, 
but unfortunately not containing the necessary material on the Fitting ideals. 



3.1 Ideals and varieties 

Let K be an algebraic field and let K[ ] be the ring of polynomials in xi, . 

with coefficients in K. A subset / C K[xi, . . . , x„] is an ideal if it satisfies 

(i) G /. 

(ii) lff,gel,thenf + gel. 

(in) If / G / and h G K[xi, . . . , then hf G /. 

Ideals are often given by generators. Let /i, . . . , /s G K[xi, . . . , Then the set 



(/i, ■■■Js) ■■= \ ^ hifi \hi,...,hse K[xi, 



i=l 



is an ideal generated by fi, . . . , fs- Any set of generators is called a basis. 

Ideals are purely algebraic objects. The geometrical counterpart of an ideal is its locus, 
or variety. Let / be an ideal in K[xi, . . . , Its corresponding variety is 

Vf(/) = {(ai, . . . , a„) e F" I /(ai, . . . , a„) = V/ G /} 

where F is some field extension of K. Note that it is often natural to choose F different 
from K. If the field is clear from context we will sometimes write simply V(/). 

Now different ideals may have the same variety. However, if one is interested mainly 
in the variety then it is useful to define 

V7 = {/ G K[xi, . . . I /" G / for some n > l}. 

If I is an ideal, then a/7 is the radical of /; it is the biggest ideal that has the same variety 
as / and all ideals having the same variety have the same radical. Also, always I C a/T 
and if J = y/l we say that / is a radical ideal. Some rudimentary properties among ideals 
and their varieties are in the following 

Lemma 3.1. Let / and J be ideals. Then 

1. v(/u J) = v(/)nv(j). 

2. v(/n J) = v(/)uv(j). 

3. / C J if and only if V(/) D V(J). 

Next we have to express the rank condition algebraically. To this end we need 

Definition 3.1. If J = (/i,...,/^), its Fitting ideal Fj is the ideal generated by all 
maximal minors of the Jacobian matrix of (/i, . . . , /s)EI 



^In general one can define Fitting ideals of minors of any given size. However, the above definition is 
sufficient for purposes of the present paper. 



Now V(F/) corresponds to the points where the rank is not maximal. However, the 
points are required also to be on V(J). Hence we conclude that the set of singular points, 
S, is given by 

s = y{iu Fi) 

In analysing varieties it is often helpful to decompose them to simpler parts. Similarly 
one may try to decompose a given ideal to simpler parts. This leads to following notions. 

Definition 3.2. A variety V is irreducible iiV = ViUV2 implies V = Vi oi V = V2. 
An ideal / is prime if f,g E K[xi, . . . , x„] and fgEl imply that either f E I oi g E I. 

There is a very close connection between prime ideals and irreducible varieties. The 
precise nature of this depends on the chosen field. However, for our purposes the following 
is sufficient. 

Lemma 3.2. If / is prime, then V(J) is irreducible. 

Any radical ideal can be written uniquely as a finite intersection of prime ideals, 

V7 = Ji n ■ • • n J^, 

where Jj ^ Ij for i ^ j. 

This is known as the prime decomposition of vT and the Jj's are called the minimal 
associated primes of /. The above Lemma then immediately gives: 

Corollary 3.1. 

V(/)=V(v^) = V(/i)U---UV(/,), 
where all V(/j) are irreducible. 

Hence our strategy in analysing varieties is to compute the minimal associated primes 
of the relevant ideal, and then examine each irreducible component separately. 

3.2 Grobner bases 

An essential thing is that all the operations above, especially finding the radical and the 
prime decomposition can be computed algorithmically using the given generators of /. To 
do this we need to compute special bases for ideals, called Grobner bases. We will only 
briefiy indicate the relevant ideas and refer to [CL092j and [GP02j for more details. 

First we need to introduce monomial orderings. All the algorithms handling the ideals 
are based on some orderings among the terms of the generators of the ideal. 

Intuitively, an ordering >- is such that given a set of monomials (e.g. terms of a 
given polynomial), >- puts them in order of importance: given any two monomials := 
x^*^ . . . x"" and x^, where ct ^ P are different multi-indices, then either x" >- or x^ )^ x". 
A common choice is to use degree reversed lexicographic ordering |CL092j . In our analysis 
we shall frequently need product orders, which are formed as follows: if >~a and >~b are 
two orderings, we shall divide the variables Xj into two subsets, and use >~a on the first 
subset and >~b on the second. This is indicated with the following notation: 

K[(X4,X5,X7), (xi,X2,X3,X6)]. 

This is the same set as K[xi, . . . ,xj] but now the parenthesis indicate that we will use 
among the variables (X4, X5, X7), and >~b among the variables (xi, X2, X3, xq), and moreover 



all monomials where variables of the first group appear are always bigger than monomials 
where there are only variables of the second group. We will see later why this is useful. 

Finally, the aforementioned Grobner basis is a special kind of generating set, with 
respect to some ordering. Given any set of generators and an ordering, the corresponding 
Grobner basis exists and can be computed. The relevant algorithm is usually called 
the Buchberger algorithm. The drawback of this algorithm is that it has a very high 
complexity in the worst case, and in practice the complexity depends quite much on the 
chosen ordering Jl 

Anyway Grobner bases have proved to be very useful in many different applications. 
Nowadays there exist many different implementations and improvements of the Buch- 
berger algorithm. We chose to use the well-known program Singular [GPS05] . |GP02] in 
all the computations in this paper. 



4 Analysing singularities 

4.1 Geometric description of the singularities 

Now getting back to our system ([3]) we see that we can take the components of p to be 
elements of Q{a,b,w)[c, s] where Q{a,b,w) is the field of rational functions of a, b, and 
w. Hence we have an ideal J = (pi,...,pi3) C Q{a,b,w)[c, s] and the corresponding 
Fitting ideal Fj. On the other hand we may view the "parameters" a, b, and w also as 
variables since they appear polynomially in the equations; hence we could also consider 
J C Q[a, b, w, c, s]. Taking this point of view we can give an intuitive description of what 
kind of situations we can expect. 

J C Q[a, b, w, c, s] 

In this way \r{J) should be 12 dimensional (recall J is generated by 13 equations), i.e. 
a curve depending on 11 parameters. On the other hand if we fix parameters a, b, and 
w we get a curve in M}'^ which will be denoted by Va,b,w In the same way we can view 
Vm(JUFj) as a variety in M^^, and fixing the parameters we get the singular points V/^^. 
Obviously V^%^ C Va,b,^ C 

Then what kind of variety should Vk(J U Fj) be? Since the Jacobian of p is of size 
13 X 14, generically we expect to get 2 independent conditions in order the rank to drop. 
That is, augmenting J with Fj should bring in 2 more equations. Hence we expect that 
Vir( J U Fj) is 10 dimensional; in other words we expect that if 11 parameters are chosen 
independently then V"/^ ^ should be empty. On the other hand if a single condition among 
parameters is satisfied, then Vfi^ ^ should consist of isolated points. 

Further, if there are 2 conditions among parameters (i.e. 9 parameters freely chosen), 
then it would be possible that V/j, ^ were one dimensional. But then our original constraint 
equations would be redundant, i.e. there would be more than one degree of freedom. 

Below we will in fact observe that if a certain condition on parameters is satisfied, 
Vfjj ^ is indeed a finite set of points. 



^So far, no satisfactory theory of Grobner basis complexity has been done. 



4.2 Singular variety 



To study V]r( J U Fj) we could in principle use Grobner basis theory in a straightforward 
manner. Let G be the Grobner basis of J U Fj using the product order Q[(c, s), (a, b, w)]. 
Let us denote hj gi, . . . , gr the elements of G which do not depend on c and s. 

Definition 4.1. Let Sj = {gi, . . . , gr); then we say that V]g(S'j) C is the singular 
variety associated to J. 

It follows from the Grobner basis theory that Va^b,w can have singularities only if 
(a, b, w) e Vk(5'j). Hence theoretically, we could now find the singularities of the Andrews' 
system in a straightforward manner by calculating the Grobner basis of JUFj. But this is 
an enormous task, due to Fj being generated by high degree polynomials, not to mention 
including the 11 parameters a,b,w. We could not get the solution in a finite time using 
our work station with 64GB memory. 

Instead, something else needs to be done. Luckily there is another approach: noting 
that pi,P3,P5 have common terms, as well as P2,P4,P6, gives us motivation to study two 
subsystems. One spanned by P5—P3 and Pq—Pa, the other one spanned by p^—pi and PQ—p2 
(along with the relevant trigonometric identities from pr, . . . ,^13)- These subsystems are 
handleable and give useful information for the whole system as well. Proceeding in this 
way we could at least determine that the singular variety is not empty and we could 
compute some subvarieties of it. 

4.3 Subsystem 4567 

Intuitively, the nodes and bars 4, 5, 6, 7 formulate a subsystem, see Figures [1] and [2l We 
suspect that when the lengths 04, . . . , 07 are such that the "4567" system is able to become 
one-dimensional, hence in some sense degenerated, there should be a singularity in the 
whole system (see also the net example in |Arp01| ). We will shortly see that this is indeed 
the case. 
Define 

<?! := P5-P3 = 04(5405 + C4S5) + asCs - aQ[cQCr - SQS7) - a^s^ 

q2 ■= P4-P6 = a4:{ciC^ - S4S5) - agSs + a6(s6C7 + CqSt) - a-jC-j 

Qi := pi+7 = c^+i + Si+i - 1, 2 = 3, . . . , 6. 

Note that qi, q2 contain only angles q, Si and parameters Oj for i = 4, . . . ,7. That is why 
we do not need the other pj's. Let J4567 be the ideal spanned by gi, . . . , gg- Hence we have 

-'4567 C Q[(C4, S4, C5, S5, Ce, Sq, C7, S7), (04, 05, Og, 07)] (4) 

where we have indicated the relevant product order. The Grobner basis G for Ji^erUFj^^^^ 
with respect to this ordering contains 191 elements (denoted by (71, ... , giQi), out of which 



3 are especially enlightening: 



916 = 040405, and 





g 


9i = 


TTtj, where 

i=l 


h 


= O4 — O5 — Og — O7 




= 04 — 05 + 06 + 07 




= 04 + 05 + 06 + 07 




= 04 + 05 — 06 — 07 


h 


= 04 — O5 + 06 — O7 




= O4 — O5 — Og + O7 




= 04 + 05 — 06 + 07 


h 


= 04 + 05 + 05-07 



Since gi is the only generator which does not contain any variables q and Si we conclude 
that 

Theorem 1. The singular variety of J4567 is 

Sj,,,,=y{{g,)). 

Note that the factorization of gi gives us the prime decomposition of {gi) and hence 
decomposition of y{{gi)) into 8 linear irreducible varieties. 

Our next task is to show that at least some points of the singular variety extend to 
actual (physically relevant) singularities of the whole system. Recall that each generator 
gi corresponds to an equation gi = 0. Since Oj > in physically relevant cases, generators 
g^ and giQ imply that all the singularities of J4567 have necessarily C6 = C4 = (conditions 
for the angles 4 and 6). In other words, in ideal-theoretic language, we can as well study 
the ideal 

T := (J4567, -^J4567' C4, Ce). 
Now the prime decomposition of ^/T has 16 components: 

= Ti n . . . n Tie. (5) 

Inspecting the generators of each of Tj, it is noticed that every Tj contains the tj's or Oj's. 
Recall that a generator Oj in an ideal corresponds in the variety to a condition Oj = which 
is non-physical. Moreover, ^3 is now a non-physical condition contradicting Oj > OVi. 
Hence we discard (as in |Arp01| ) those ideals which have a non-physical generator that 
would imply Oj < for some i, and we are left with 7 ideals, whose generators are: 





= (c? + - 


1, 


tl, 


se + 1, 


■55 ~ C7, 


C5 + S7, 


S4 + 1, C4, Cq) 




= (c? + - 


1, 




se + 1, 


■55 + C7, 


C5 - ^57, 


S4 + 1, C4, ce) 




= (c? + - 


1, 




■56 + 1, 


■55 + C7, 


C5 — ■57) 


S4 — 1, C4, Ce) 


T4 


= (c? + - 


1, 




Se - 1, 


■55 - C7, 


C5 + ■57, 


S4 + 1, C4, Ce) 


T5 


= (c? + - 


1, 




Se - 1, 


■55 + C7, C5 - S7, 


S4 + 1, C4, Ce) 




= (c? + - 


1, 




■56 - 1, 


■55 — C7, 


C5 + S7, 


S4 — 1, C4, Ce) 


T7 


= (c? + - 


1, 




■56 - 1, 


■55 + C7, 


C5 - ■57, 


S4 - 1, C4, Ce) 



Especially, we see that Sq = ±1, S5 = ±07, C5 = ±57, and S4 = ±1. Now we are ready to 
continue with the original system J U Fj. 

Remark 4.1. Mathematically speaking the analyses of all cases Tj are completely similar. 
However, on physical grounds the cases Ti, T2, Tg and T7 are not so interesting. Indeed, in 
these cases the length of one of the rods corresponding to 04, 05, ag and aj is equal to the 
sum of the lengths of three others. Hence all four rods could be modelled as a single rod 
which would make the whole model significantly simpler. In the remaining cases no such 
reduction can be done, and we chose to examine the ideal T5 in detail. See also remark 

m 

The case T5 gives us conditions S4 = —1, sq = 1, S5 = — C7, C5 = S7, and ay = a^+aQ—a^ 
which we substitute into the original system. Next we will show that the resulting system 
has real solutions. These will be the required singular points. 

The above substitutions simplify the generators of JU Fj so that we get the following 
ideal: 

K= {K,UK2), 

' ki = a2(-CiC2 + S1S2) + Citti - 3303 - bi 

a2(-SiC2 - C1S2) + Sitti + Cstts - 62 

1 

1, (6) 
57(04 — as) + 5303 + bi — wi 
07(05 - 04) - csas + b2-W2 



k2 

ki 

(k, 

ka 
kr 
h 



cl + sl 

^2 + ^2 



4 + sl-i 



c^ + s^- 1. 

In K2 we have 4 equations for 4 unknowns C3, S3, Cj, and Sj] hence it appears reasonable 
that we can get a finite number of solutions. Then we can substitute the computed values 
to Ki which then becomes also a system of 4 equations for 4 unknowns ci, si, C2, and 
S2- By the same reasoning we again expect that it is possible to get some solutions for 
appropriate parameter values. 

We could numerically solve the variables from these equations (and, indeed, we will, 
in the numerical examples), but to analyze the situation in more detail we need to study 
these further. 

Then starting with the system K2 we solve the angles 3 and 7 by the following trick. 
First we inspect the ideal generated by K2 in the ring 

Q(6i, 62, wi, W2, as, 04, as) [03, S3, C7, sy]. 

Calculating the Grobner basis G of {K2) with respect to the lexicographic ordering we get 
4 generators: 



91 = /1S7 + /2S7 - hfi 

92 = 2(62 - w^2)(a4 - 0-5)07 - 2(61 - wi){a4 - as)s7 + /s 
^3 = 0.383 + (04 - 05)37 + bi-wi = 

h = 03C3 + (04 - 05)07 + W2 - 62 = 0. 



(7) 



where the auxiliary expressions /j are lengthy combinations of the parameters a^, bi (see 
the appendix) I2I 



•^The algorithms actually give by default only sums of monomials instead of products like 2(62 — 



Now gi contains only S7 and parameters. Note that /i = if and only if 
Assuming 04 7^ 05 the equation ^1 = is a polynomial in S7 of degree 2, hence in order to 
have real solutions we need to impose the condition 

/2' + 4/1/3/4 >0. (8) 

This condition can easily be checked when the parameters a, b, w have been given numerical 
values. Once S7 is known, C7, S3, C3 can be solved from the linear equations of G, provided 
04 7^ 05 and W2 7^ &2- 

The cases W2 = &2 and/or = can be summarized as follows: 

(i) If W2 = &2 but a4 7^ as, we still get equations similar to G, but now S3 has a quadratic 
equation instead of S7. 

fii) If (34 — 05, the system typically does not have solutions. At least, a further condition 
among parameters, namely |6 — w| = 03, arises. We shall not elaborate this non- 
generic behaviour further. In Section l4.5.2l we consider an example of this situation. 

Remark 4.2. In general, when the inequality in ([8]) is strict, S7 has 2 possible values. 
Therefore, the tuples (S3, C3, S7, C7) have in general 2 possible values because the other 
ones in the tuple are determined uniquely from S7. 

The only thing left to be done, in this J4567 subsystem case, is to solve Ci, Si, C2, S2. This 
is done with the ideal (Ki) given in ([6]). 

Remark 4.3. Had we used any other Tj instead of T5 above, we would have ended up with 
this same ideal (Ki). 

We calculate the Grobner basis G of {Ki) , this time in the ring 

Q(ai, 02, 03, bi, &2, C3, S3) [ci, si, C2, S2]. 

Note especially that S3, C3 are here treated as parameters, due to being now known expres- 
sions in the parameters a, b, w. We again use lexicographic ordering and get 4 generators 
gi,...,g4. Analogously to S7 above, now for S2 we get the second degree polynomial 
equation 

gi = (— 4a]^a2)s2 — nin2 = (9) 

where 

rii = al + 2aia2 + — al — 2a3biS3 + 2036203 — b\ — b\ 
n2 = a\ — 2aia2 + — — 2a36iS3 + 2036203 — b\ — b\ 

and linear equations for C2, Si, Ci: 

g2 = diC2 + d2 + d3 
h = ksi + l2 + h 

where the auxiliary expressions di, U are certain known (but lengthy) functions of a, 6, 
apart from l^^ which depends on Si,S2,C2 as well. (See the appendix.) In order to have 
real solutions for S2, ([9]) imphes the condition 

E := nin2 < 0. (10) 

W2){0'A — as) but we have simplified these by hand. Also Singular [GPS05| could be used to automatically 
factorize into products but would involve some more elaborate programming. 



These Qi determine S2,C2,Si,Ci provided di ^ 0, li ^ 0, ai 7^ 02- To analyse the cases 
di = 0, Oi = a2, and/or /i = 0, it is helpful to define 

do := al + 2036153 - 2036203 + 6^ + 62. 

It turns out that li = 4^ di = <^ do = 0. After rearranging the terms (see the 
appendix) it can be seen that the condition (fTOl) is equivalent to 

(ai - 02)^ < c^o < (oi + 02)^- 

Therefore, if Oi 7^ 02 then lio 7^ and the equations above can be solved. The case ai = 02, 
do ^ does not essentially change the situation: we still have a quadratic equation for S2, 
and linear ones for the others, with a different coefficient for ci. 

The remaining case ai = 02, c?o = corresponds to the situation where the centre node 
coincides with the origin. This gives another singularity (the angle yi remains arbitrary) 
but is a rather special case and will not be pursued further here. 

Theorem 2. Let us suppose that the parameters a, 6, w satisfy the following conditions: 

04 7^ as and 

ni{Aaia2 - rii) > 

f^ + lQ{a,-a,f\b-w\'hh>0 & 

Then Va^b,w contains at least 2 singular points. If the inequalities are strict we get in 
general at least 4 singular points. 

It may appear that we also have at most 4 singular points. However, it is a priori 
possible that the other systems Tj yield more singular points with the same parameter 
values. 

Proof. The first part of the theorem merely collects what we have shown above, with the 
simplifications n2 = rii — 4aia2 and /i = 4(04 — a^Ylb — wp. The conditions are due to 
univariate second degree polynomial equations, which have real solutions if and only if 
([8]) and (fTOj) (for sj and S2, respectively) are fulfilled. The other variables are determined 
from linear equations: S4, C4, . . . , Sq, Cq from T5; S3, C3, C7 from Ki, Si, Ci, C2 from K2. 

For the number of singular configurations, note that we have second order equations 
for S7, hence at most 2 values for the tuple (53,03, 57,07), and 52. So in general if there 
are two separate roots both for 57 and S2, we get four different singularities. □ 

Similar results can be presented for any Tj but we will not catalogue them here. 
4.4 Subsystem 367 

Comparing to examples in [ArpOl] it was perhaps intuitively clear that subsystem J4567 
produces singularities. It is a bit more surprising that there is another subsystem produc- 
ing singularities: the one formed by the nodes 3, 6, and 7. 
Define 

hi ■■= -P5 +Pi=ae {cecj - sqSt) + a^s^ - a^s^ + wi - bi 
h2 ■= -P6 +P2 = a(i [secj + cqSj) - arCj + 0303 + W2 - 62 
h3 := P9 = C3 + S3 - 1 
hi := P12 = Cg + Sg - 1 
h5 ■= Pl3 = C7 + S7 - 1. 



It is important to note that hi, h2 contain only angles 3,6, and 7, therefore only pg, pi2, Pis 
are relevant to them. As parameters we now have not only the lengths 03, Og, 07, but also 
bi, . . . ,W2 i.e. the positions of the fixed nodes A and B in Figure [2l Let J^qj be the ideal 
generated by /ii, . . . , h^. We will proceed in a similar way as with the subsystem J4567. 

First we will consider the singularities of the subsystem J367 using the following product 
order: 

-^367 U Fj3g7 C Q[(c3, S3, Cg, se, C7, S7), (03, ag, ay, 61, 62, ^i, W2)] (11) 

The relevant Grobner basis G contains 96 generators of which two are especially interest- 
ing: 

4 

gi = Y\ Zi where 

i=l 

Zi = (cts — cte + (^lY — \b — (l^) 
Z2 = (0.3 + 06 + ajY -\b - 
^3 = (0-3 + - a-rY - |6 - 
Zi = (03 — — a-jY — \b — wf'. 
The latter one gives us the singular variety Sj.^^^. 

Theorem 3. The singular variety of J367 is 

Remark 4.4. It is worth noting that, contrary to the linear constraints in Theorem [1] 
related to J4567, the Zi in Theorem [3] give quadratic constraints = related to J367 
and have the interpretation "|a3 ± ag ± 07] = distance between the fixed points A and 
B". Furthermore, again the factors Zi give the irreducible decomposition of the singular 
variety. 

Since > 0, we get Cg = from gi2 = 0. This simplifies computations considerably. 
Let us define 

U := (^367, -^J3675 Cg). 

The prime decomposition of U turns out to have 8 components: 

Vu = Uin---nUs. 

Inspecting the generators of each of Ui, it is noticed that the ideals Uk, k = 5 . . .8 contain 
generators which imply = for some i. Hence those are discarded as non-physical and 
we are left with 4 ideals: 

Ui = {Ui, U2, + Sj- 1, Cg, Sq - 1, S3 + S7, C3 + C7) 

U2 = {Ui, U2, Cj + Sj- 1, Cg, Sg + 1, S3 + S7, C3 + C7) 

Us = {Ui, U2, + Sj- 1, Cg, Sg + 1, S3 - S7, C3 - C7) 

f/4 = U2, C7 + S7 - 1, Cg, Sg - 1, S3 - S7, C3 - C7) 



ui = -SQCjaQ - C3a3 + C7a7 + 62-^^2 
U2 = SQSjaQ + S3a3 — S7a7 + 61 — Wi. 



With these, we continue studying the whole system JUFj. Each Ui will lead to a different 
case with sq = ±1, S3 = ±57, C3 = ±07. Let us look for example the ideal t/iQ This gives 



S6 
C7 

C3 
S3 



62 - W2 
Clg — CI3 — CI7 

bi - Wi 

-C7, 
-S7. 



(13) 



We should expect to run into an equation Zi = for some i, where the expressions Zi are 
given in ( fT2l) . Combined with C7 + — 1 = the equations (fT3l) give = 0. Likewise, t/j 
implies Zi = for i = 2, 3, 4. 

Remark 4.5. The condition 2:2 = is physically a redundant case: it means that the 
system can barely reach from A to B when the subsystem of the rods 03, ag, 07 is fully 
stretched, i.e. it has no room to move. Therefore also U2 corresponds to a rather trivial 
case. See also Remark [4. II 

Using Ui we can now eliminate the variables corresponding to angles 3, 6, and 7. Doing 
the substitutions in J U Fj we are left with the following generators. 



L 



{L1UL2), 
' h = 02 
h 
h 

M 

(h 

k 
h 



C1C2 + S1S2) + Cifli + ^703 

a2(— S1C2 — C1S2) + sitti — cja^ 

Co + Sn 1, 



hi 
h2 



(14) 



04(5405 + C4S5) + Csag + 87(06 - 07) 

04(0405 - S4S5) - 5505 + 07(00 - 07) 



C4 ~l~ Sa 



ci + si 



where the S7, C7 are no longer variables, but known expressions from ( fT3l) and kept here 
only for clarity of notation. 

Remark 4.6. Before working on Li and L2 we comment briefly on the other Ui cases. 
Introduce L3 and L4: 



Lq : < 



02(— C1C2 + S1S2) + CiOi — S7O3 
02(-SlC2 - C1S2) + SiOi + C7O3 

C2 + - 1 = 



&1 
&2 



La : < 



04(5405 + C4S5) + C5O5 - S7(06 + O7) 
04(0405 - S4S5) - S5O5 - C7(06 + C7) 

cl + sl-l 
cl + sl-l 














0. 



*As with J4567 and T5, the other cases are completely similar and we will comment them shortly. 



Had we used U2 instead of Ui, we would end up with the system Li, L4. Likewise, U3 
would give the system L3, L2, and U4 would give the system L3, L4. Yet another point of 
view is, that Sq = ±1 picks between L2 and L4, while (03,53) = ±{cj,Sj) picks between 
Li and L3. More precisely, = 1 {sq = —1) gives L2 {L4), and (03, S3) = (— C7, —sj) gives 
Li. The choice (03, S3) = (07, S7) would give L3. 

Continuing with Li and L2, we notice that L2 contains only the variables C5, S5, C4, S4 
(angles 4 and 5), has 4 equations and 4 variables hence is expected to have a finite solution 
set and will be handled analogously to the ideal K2 in ([6]). Calculating its Grobner basis 
G in the ring 

Q(a4, as, ae, 07) [(04, C5, S5, C7, S7), (54)] 
we obtain 12 generators, the first one being 

gi = 2040584 + 04 + Og — Og + 20507 — O7. 

Hence S4 can be explicitly solved: 

o| + O5 — Og + 2O6O7 — O7 

S4 = . 

— 204O5 

The other generators are too messy to be of much use. Then using the formula cl 
we get 



(15) 
= l-sl 



(04 + O5 — Og + 07)(04 — O5 + Og — 07)(04 — O5 — Og + 07)(o4 + O5 + Og — O7) 

4o|o| 

— XX- (16) 



40405 

The product term in the numerator has to be nonpositive, in order to have any real 
solutions: 

hW^ < 0. (17) 

After solving S4, C4 we can proceed to solve S5 and C5. For this we use the ordering 

Q(04, O5, Og, 07)[C5, S5, C4, S4, C7, Si] 

and pick the two relevant equations from the corresponding Grobner basis: 

(— Og + 07)85 — a^c^Sj + ttiSiCi + 05C7 = 
( — Og -|- 07)05 — O4C4C7 — O4S4S7 — ci^s^ = 0, 

which are linear equations for 85,05, provided Og 7^ 07. 

Remark 4.7. In the case Og = 07 the situation is different: L2 then decomposes into 

3 prime ideals, of which only one is physically feasible and gives a singularity only if 

04 = 05. Thence this is a rather special case and will not be considered further here. 

The subsystem L2 is now fully solved. Moving on to Li, we will see that the analysis 
is very similar to that of Ki from ([6]). Therefore we will skip some details. After forming 
the Grobner basis of Li in the ring 

, 62, fll, 02, O3, C7, S7)[ci, Si, C2, S2] 



with respect to the lexicographic ordering, we get for S2, after simphfications, the relation 



4 




(18) 



Again for the real solutions the numerator has to be nonnegative 



n3(4aia2 - ^3) > 



(19) 



We can now solve C2, si and ci, provided their coefficients are nonzero, from the linear 
equations 



and Ti are lengthy, yet polynomial, expressions in the parameters, apart from r4 which 
depends on Si, S2, C2 as well. (See the appendix.) 

What about the cases 71,4 = and/or ai = 02? It can be shown, as with d^, that the 
condition ^3(40102 — n^) > is equivalent to 



Therefore, if ai 7^ a2 then 77,4 7^ and the equations above are sufficient. The case Oi = 02, 
n4 7^ does not essentially change the situation: we still have a quadratic equation for 
S2, and linear ones for the others, with a different coefficient for ci. 

The remaining case ai = 02, n4 = is analogous to the n2 = case within J4567 and 
likewise will not be pursued further. 

Theorem 4. Let us suppose that the parameters a,b,w satisfy the following conditions: 



Then contains at least 2 singular points. If the inequalities are strict we get in general 
at least 4 singular points. 

Similar results can be represented for any V(f/() but we will not catalogue them here. 

Proof. The last two conditions are due to univariate second degree polynomial equations, 
which have real solutions if and only if fl20|) (for S2) and fl2T|) (for C4) are fulfilled. The first 
condition is needed for the other variables to be determined uniquely: S3, C3, sq, cq, sj, cj 
from V(f/i), S4, S5, C5 from L2, and Si, Ci, C2 from Li. 

For the number of singular configurations, note that we have second order equations, 
hence at most 2 values, for C4 and S2. So in general if there are two separate roots both 
for C4 and S2, we get four different singularities. □ 



2aia2n4C2 — Aa\a\s\ + ri = 0, 
- 2ain4Si + r2 + r3 = 0, 
{a\ — a\)ci + r4 = 0. 



where 



= \h\^ + a\ + 2a3(62C7 - his-j) 



(ai - 02)^ < ^4 < (ai + 02)^. 



n3(4aia2 - ris) > 

tjt^tQtg < 0. 



(20) 
(21) 



4.5 Two special cases with symmetry 

Let us look more closely at two special cases: 04 = ag, = oy, and either 04 = 05 or 
04 7^ 05- 

4.5.1 The case 7^ 05 

Motivated by the original benchmark values |Sch90] we give the following 

Lemma 4.1. When CI4 — CLq and 05 = 07, there is a relation between the angles 4 and 6: 
either i/q = —1/4 or = 2/4 + tt. Furthermore, if also 04 7^ 05, the angle yj variables, i.e. 
C7, S7, are uniquely determined from 04,^4,05,55. 

Proof. Looking for relations between solely angles 4 and 6, we substitute = og and 
as = 07 to the subsystem J45g7 and formulate a suitable elimination ideal. In ideal- 
theoretic language, we define 

ri := 04(3405 + O4S5) + 0505 - a4(og07 - 3^37) - 0557 
r2 := 04(0405 - S4S5) - 0555 + a4(sg07 + OgS7) - 0507 
= + sf^3 — 1, i = 1, . . . , 4, 

where = Qi with substitutions CI4 OiQ and 05 = 07, and investigate the ideal / := 

(ri, . . . , rg) in the ring 

Q(a4, as, ag, a7)[(o5, 35, 07, S7), (04, S4, og, sg)]. 
Calculating the elimination ideal /4,g := J fl Q[o4, S4, og, sg] we get 

hfi = (S4 + Sg, Og + Sg - 1, O4 + S4 - 1). 

Calculating the prime decomposition of ^/hfi we get 

^/hfi = (Og + Sg - 1, O4 - Og, S4 + Sg) n (Og + Sg - 1, O4 + Og, S4 + Sg). 

Since /4,g C/cJcJUFj, we have 

V(/4,g) D V(JUFj). 

From these prime ideals we can see that everywhere in V(/4_g), and therefore in the variety 
of the singularities of the whole system as well, Sg = — S4 and either Og — O4 or Og — — O4. 
These translate into two possible relations between the angles 7/4 and ?/g. 

(C6, Sg) = (04, -S4) yg = -2/4, (Og, Sg) = (-O4, -S4) ^ yg = ?/4 + TT. (22) 

This proves the first claim. If we take into account either one of the prime ideals of ^/l^fi 
in I and calculate the Grobner bases we get ideals where 07 and S7 depend linearly on 04, 
S4, os and ss, and can be explicitely solved, as we will show next to prove the latter claim 
of the lemma. For the case (sg, og) = (— S4, — 04) we get 

which imply yj = y^ + ^. (23) 
S7 = 05 2 



For the case (sg, Cg) 
cated: 



-S4, C4) the expressions are, albeit hnear, shghtly more comph- 



€7(04(54 - C4) - 05(20454 + 05)) + 57(204(05 + O4C4S4)) - S5((04 + O5) - 2O4O5S4) = 

-07(2040454) + 57(04(04 - S4) + 05)) + (04 - 05)05 - 204O5S5O4 = 0. 

We prove that these indeed determine 07, S7: all we need to do is check that the determinant 
of the coefficient matrix A of the linear equations does not equal zero: 



A :-- 



04(54 



05(20454 + 05) 204(05 + 0404^4^ 



-204O4S4 



O4 yc^ 



sf) + ai 



prove det(yl) 7^ 0. 



Now det{A) simplifies due to O4 + S4 = 1, resulting in 



det{A) = 20405(04 + 05)(04 — 05)54 + (04 — 05)(04 + 05)(04 + O5) 

Let us then consider det{A) as a function of S4. Since S4 G [—1, 1], det{A) : [—1,1] 
Clearly if 04 = 05, det{A) = so we need to assume 04 7^ 05. Set 



h{s4) 



det{A) 



(O4 + 05)(04 - O5) 



204O5S4 + (04 + O5) 



and inspect when h = 0. Since 04 > and 05 > the linear function h has its minimum 
at -1. 

h(—l) = O4 + O5 — 204O5 = (05 — 04)^ > 0. 

This proves h always, therefore under the assumption 04 7^ 05 also det{A) 7^ as 
claimed. □ 



4.5.2 The case 04 



05 



We study the special case 04 = 05 = 05 = 07, whence the 4567-subsystem is capable of 
"buckling" in more complicated ways, thereby producing further interesting configurations. 
This resembles then the net example in |Arp01| . 

Let us see how J4567 simplifies with substitutions 04 = 05 = 06 = 07. Note that the 
assumptions of Lemma 14.11 considering t/j are no longer valid. Let 

I '■= <^4567 with 04 = 05 = 06 = 07 and S6 = —S4 
and compute its prime decomposition. This results in 

Vl = Ji n /2 n /a with generators 



sl + cj- 1, 

04 — 06, 

0? + - 1, 
S5 + O7S4 

05 — O7O6 



S7C6, 
S7S4 



C65 

S4 + 1, 
O4, 

Cy -\- Sj 

4 + 4 



4 + 4 

C4 + Cg, 

Cj -\- 



S5 
C5 



C7, 
S7 



(24) 



Each of these has a geometrical interpretation, see Figure [3l I2 corresponds to 1/4 = 
— 7r/2,|/6 = 7r/2 which means that nodes A and P2 coincide. This is like the T5 situation. 
Indeed, the ideal J U Fj U I2 turns out to be exactly T5 with the extra condition 04 = 05. 



Although it is not immediately apparent but in that situation there also arises a new 
condition among the parameters: 03 = |6 — i.e. "03 equals the distance between A and 
B". Note that here the Fitting ideal Fj^^^^ has not been used at all, contrary to the T5 
calculations. 

I3 corresponds to ?/6 = ?/4 + tt and 1/5 = 7/7 — 7r/2 so that now nodes P3 and P4 coincide. 
Then again, Ji corresponds to i/q = —y^ and 1/5 = ye + yj, which interestingly is not a 
singularity but merely expressing a symmetry in the system due to 
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Figure 3: The configurations corresponding to Ji, J2, /3 in the case 04 = 05 = oe = 07. 
4.6 Other subsystems 

Now contemplating Figure [2] we see that it would be possible to find other singularities by 
analysing still other subsystems. For example the subsystem corresponding to rods 3, 4 
and 5 is by symmetry similar to subsystem 367: we simply exchange the roles of variables 
and parameters associated to rods 4 and 6, and 5 and 7. Further we could consider other 
subsystems formed from different "paths" between the nodes A, i?, O: i.e. subsystems 
</i23, -^^1245, </i267- Again by symmetry the system J1267 is completely similar to J1245, but 
cases J123 and J1245 give new singularities. We checked that in these cases the singular 
variety is not empty, and that at least for some parameter values we get singular points. 

We did not analyse these cases in detail because computations are quite similar to 
those given above for subsystems J4567 and J367. Hence we did not feel including these 
would give significant additional value and therefore left them out to avoid expanding this 
quite a long presentation further. 

5 Numerical examples 

In this section we will calculate numerical examples for both types of singularities. In- 
terestingly, the explicit expressions within G, G, as well as in the Grobner bases of Li 
and L2, are unstable for numerical computations. It is better to use the original defining 
equations of -f^i, i^'2; -^i? -^2 in the computations. We shall not explore this stability issue 
here due to its non-relevance for the present context. 
We present 4 examples: 

1. The original benchmark parameter values, see |MI03] . We show that then the system 
is avoiding singularities^ 

■^Thereby validating its benchmark status. That is, the numerical difficulties encountered there are 
indeed due to the "numerical stiffness" of the problem, not to a nearby singularity. 



2. We explore how should ai, 02 be changed in order to have J4567 type singularities in 
the system. Here we have an interpretation for the result: the lengths Oi, 02 must be 
such that the subsystem 4567 has room for a certain kind of "buckled" configuration. 

3. We explore how should 61, ai, 02 be changed in order to have J367 type singularities 
in the system. 

4. A special case which shows a rational solution, that is Cj, Sj e Q for all i. This shows 
unambiguously that we can find singular points because in this case there are no 
numerical errors related to floating point computations. 

5.1 Original values 

In this example, we will use the original values for the parameters Oj, 6j and show that the 
system then has no singularities. The original parameters used in the benchmark tests 



jSaiQOl IHW911 [Ml03] are 

ai = 0.007 02 = 0.028 ag = 0.035 = 0.020 05 = 0.040 ag = 0.020 = 0.040 



Since ay = 05 and Og = 04, we have ^4 = tg = (and ti < 0, < 0) so we could have an 
J4567 singularity: T3 or T5. 

Remark 5.1. Interpretation: both T3 and T5 describe a situation where the 4567 system 
has 'collapsed' into a 1-dimensional object. The ideal K2 tells us how restricts the 
possible attitudes of 4567. In T5 the centre node P2 has been pushed in, in T3 it has been 
pulled out. 

Let us look more closely first at T5, say, and check the conditions (E]) and (ITOl) . The 
first one is fulfilled. For E we first need to solve C3, S3 from \J{K2). Their solutions are 

(c3, S3, C7, S7) e {(0.4299535996, -0.9028509856, -0.9975812008, 0.06951077517), 

(0.9266735994, -0.3758670513, -0.1283212011, 0.9917326602)} (26) 

With these C3, S3 we can compute E. Both sets in fl26|) give E = O{10^^) > and the 
condition (fTOj) is violated, hence there are no ( J4567— )singularities. What about other 
singularities? This is answered by the following 

Theorem 5. With the original benchmark parameter values (125|) . the Andrews' squeezing 
system has no singularities. 

Proof. We now have = Og, = aj and 04 7^ 05. Lemma WA\ implies variables Cg, Sg, 
C7, S7, and so and 1/7 can be explicitely solved in terms of C4, S4, C5, and S5. It is then 
possible to reduce the original system of constraint equations, by forgetting the last two 
equations from ([2]), and consider 

' ai cos(?/i) - 02 cos(?/i + ^2) - 03 sm{y^) - hi =0 



h 



'1 = 



-0.03635 62 = 0.03273 Wi = -0.06934 W2 = -0.00227. 



(25) 



ai sin(yi) - 02 sin(?/i + 1/2) + 03 cos(?/3) - 62 
ai cos(?/i) - 02 cos(yi + 7/2) - 04 sin(y4 + 1/5) - cos(t/5) - Wi 
^ai sin(?/i) - 02 sin(?/i + 7/2) + ^4 cos(y4 + 2/5) - ag sin(?/5) - W2 











0. 



These are equivalent to 

' ai cos(?/i) - fla cos(?/i + 7/2) - % sm(y3) - bi =0 

ai sin(yi) - a2 sin(yi + ^2) + 0,3 003(1/3) - &2 =0 

-04 sin(y4 + 7/5) - as cos(?/5) + 03 sin(?/3) + (61 - wi) =0 

^ 04 cos(2/4 + 2/5) - 05 sin(?/5) - 03 cos(?/3) + (62 - ^^2) =0 

These can be again represented as polynomials. 

777,1 := CLiCi - a2{ciC2 - S1S2) - 03*3 - h = 

m2 := fliSi - a2(siC2 + C1S2) + 0303 - 62 = 

7713 := aiCi - a2(ciC2 - S1S2) - a4('54C5 + C4S5) - asCs - Wi = 

m4 := aiSi - 02(5102 + C1S2) + 04(0405 - S4S5) - 0555 - W2 = 

mi+4 := + Si - 1 = 0, z = 1, . . . , 5. 

Substituting the original parameter values (|25l) . as rational numbers, into the polynomials 
rrii we form an ideal I := (mi, . . . , mg). Let i^' := / U Fj, where Fj is the Fitting ideal of 
J, and inspect K in the ring 

Q[(Cl, Si, C2, S2), (C3, S3, C4, S4, C5, S5)]. 

Now it is possible to compute the Grobner basis Gk for K explicitly (unlike for J U Fj in 
the introduction) and results in 

Gk = (1). 

This implies V(K) = 0, proving that with these original parameter values there are no 
singularities. □ 

5.2 J4567 singularity: original values, apart from ai,a2 

Let us see how changing ai and/or 02 might produce J4567 type singularities. Our analysis 
reveals that by suitable combinations of Oi and 02 we can get between zero and four 
singularities (of type J4567, that is). The number of singularities is determined by 03,53, 
and E. 

Considering E as a. function of ai,a2 we plot the area where E < 0. Recall that 
E depends on C3 as well, and C3 has two possible values so we get two functions: E = 
Ei{ai, 02) (resp. E = E2{ai, 02)) corresponding to the first (resp. second) value of C3 from 
(l26l) . See Figure H] where the areas inside the rectangular areas are Ei < 0. 

• no singularities: Ei > 0, £'2 > 0. 

• 1 singularity: Ei = 0, i?2 = 0, which leads (with T5) to two possible values: 

(ai = 0.05986, 02 = 0.01035), (ai = 0.01035, 02 = 0.05986) 

• 2 singularities: one of -Ei, £'2 is < 0, the other one > 0. 

• 3 singularities: one of Ei,E2 is < 0, the other one = 0. 

• 4 singularities: Ei < 0, E2 < 0. 




al al 

Figure 4: The rectangular lines are Ei = (thick line) and E2 = (thin line), the areas 
inside each Ei = line are where Ei < 0. Left panel: T5 case, right panel: T3 case. 



For example, let us concentrate on T5 and choose ai = 0.03, 02 = 0.055, say, whence 
the system is able to reach four singular configurations (see the left panel of Figure H]). 
Now Si, Ci for z = 1, 2, 3, 7 are determined by V(i^). The other values, for angles 4,5,6, are 
determined by ^{T^). The results are in the Table [H The corresponding configurations 
are visualized in Figure [H Doing similar tests with T3 instead of T5 yields the Ei areas 
in the right hand panel of Figure |H Singular configurations implied by T3, with choices 
Oi = 0.06, 02 = 0.06 which imply 4 singularities, are in Figure [61 To save space we have 
not tabulated the actual values of the angles in T3 case. 



5.3 J367 singularity: original values, apart from 6i,ai,a2 

A necessary condition to have a J367 type singularity is at least one of the Zj's vanishes 
(fT2l) . Substituting the original parameter values we notice that none of these is zero. Let 
us then investigate how we should change some of the parameters in order to have Jsr? 
type singularities. Take bi and Ui, say, and choose bi := —0.026913593 so that zi = 0. □ 
We seek to further fulfil the sufficient requirements by Ui. 

n3(4aia2 -ns) > (EO]) 
t7t5t6ts < 0, m 

and use Li,L2 to find the actual singular configurations. With the original parameter 
values te = 0, therefore ( l2TI) is fulfilled. Therefore we only need to study ( l20l) . For that, 
we proceed analogously to Example 15. 2t treat the expression 723(40102 — n^) as a function 
of ai, 02. For that, we first need 07,57. Them we get from ( fT3ll 

62 — Wo 

C7 = — = -0.6364 

Og — 03 — CLj 

.7 = -A^^ = 0.7714. 
03 + 07 — ag 

The region of ai, 02 plane where n3(4aia2 —ns) > is shown in Figure [71 We pick a value 



^This corresponds to moving B slightly to left. 



variable 


singularity 1 


singularity 2 


singularity 3 


singularity 4 


Cl 


-0.8322 


-0.4564 


-0.1157 


-0.1038 


Sl 


-0.5544 


0.8898 


-0.9933 


0.9946 


C2 


-0.3045 


-0.3045 


0.4467 


0.4467 


S2 


0.9525 


-0.9525 


0.8947 


-0.8947 


C3 


0.4300 


0.4300 


0.9267 


0.9267 


S3 


-0.9029 


-0.9029 


-0.3759 


-0.3759 


C4 














S4 


-1 


-1 


-1 


-1 


C5 


0.0695 


0.0695 


0.9917 


0.9917 


S5 


0.9976 


0.9976 


0.1283 


0.1283 


C6 














S6 


1 


1 


1 


1 


C7 


-0.9976 


-0.9976 


-0.1283 


-0.1283 


S7 


0.0695 


0.0695 


0.9917 


0.9917 



Calculating the corresponding angles we get the following values. 



Angle 


singularity 1 


singularity 2 


singularity 3 


singularity 4 


yi 


-2.5539 


2.0448 


-1.6867 


1.6747 


2/2 


1.8802 


-1.8802 


1.1077 


-1.1077 


2/3 


-1.1264 


-1.1264 


-0.3853 


-0.3853 


2/4 


-1.5708 


-1.5708 


-1.5708 


-1.5708 


2/5 


1.5012 


1.5012 


0.1287 


0.1287 


2/6 


1.5708 


1.5708 


1.5708 


1.5708 


2/7 


3.0720 


3.0720 


1.6995 


1.6995 



Table 1: The singularities of J4567 type, original values apart from ai,a2. The values are 
presented only with 4 decimals but were computed with 16 decimals. 

inside the "allowed" annulus, say ai = 0.02 and 02 = 0.055 in order to get singularities. 
Then let us find the actual singular configurations: since tg = 0, from (fT6l) we get C4 = 
and from (fT5l) S4 = —1. The other angles are found as follows: 3 and 6 from (fT3l) and 
the remaining ones 1,2,5 from L. The results are in Table [2l The corresponding singular 
configurations are drawn in Figure [HI Note that there are only two singular configurations, 
instead of four, since f|T6l) has only one (double) root C4 = instead of two separate roots. 



5.4 A rational case 

Finally, let us show a rational valued singularity, that is q, Sj G Q. Choose 

a4 = 05 = ag = 07 = 3/20 ai = 1/10 02 = as = 1/2 
bi = -1/10 62 = 1/5 wi = -2/5 W2 = -1/5 



-0.02' 

-0.08 




-0.06 -0.04 -0.02 



0.02 0.04 




-0.08 -0.06 -0.04 -0.02 



0.02 0.04 





-0.04 -0.02 



-0.08 -0.06 -0.04 -0.02 



0.02 0.04 



Figure 5: Singular positions (according to J4567, Z5) when ai = 0.03, a2 = 0.055 and 
03, . . . , 07 have the original values. One can see a physical explanation to the singularity: 
the centre node P2 is 'pushed in' so that nodes P3 and P4 coincide. 



and solve c, s from the generators of /2 U J U Fj in (12^ . Now C5, S5, C7, S7 are arbitrary 
(apart from cl + si = 1, Cj + Sj = 1) and the chosen result is (see also Figure [9]) 

c= (0,3/5,4/5,0,3/5,0,4/5) 

s = (1,-4/5,-3/5,-1,4/5,1,3/5). 



6 Conclusion 

We have studied singularities of the multibody system "Andrews' squeezing system" which 
is a well-known benchmark problem both for multibody solvers and differential-algebraic 
equation solvers. Using our tools we have shown in Theorem [5] that the original benchmark 
problem is indeed void of singularities, thereby assuring that whatever numerical problems 
in the benchmark tests are met, they are indeed due to something else than a nearby 
singularity of the system. Apparently, this non-singularity of the problem has not been 
rigorously proven in the literature. 

However, we have shown that with suitably chosen parameters (a, b, w) , this system can 
exhibit singular configurations. In fact, there are families of values (a, b, w) that produce 
singularities, see Theorems [2] and HI We provide examples of singularities, calculated 
using the original benchmark parameter values apart from bi,ai,a2- Considering ai, 02 as 
freely chosen parameters. Figures H] and [7] show the areas of ai, 02 plane where the system 



Andrews' squeezing sysli 




Andrews' squeezing syster 




-0.07 -0.06 -0.05 -0.04 -0.03 -0.02 -0.01 0.01 0.02 0.03 

Andrews' squeezing system 




Figure 6: Singular positions (according to J4567, ^3) when oi = 0.06, 02 = 0.06 and 
03, . . . , ay have the original values. One can see a physical explanation to the singularity: 
the centre node P2 is now 'pulled out' so that nodes P3 and P4 coincide. 



exhibits singularities. For example, choosing the point (01,02) within the intersection of 
the three areas in Figures H] (both panels) and [7] would give a system with 10 singular 
configurations. 

A natural question that remains is, if these presented singularities are the only possible 
ones? In other words are there singularities which do not come from the singularities of 
some subsystem? While the Grobner bases techniques in principle provide a way to answer 
this question directly, we could not do so in practice due to complexity problems. 



0.05 



0.04 - 



a2 0.03 



0.02 



0.01 - 




0.01 0.02 0.03 0.04 0.05 

al 



Figure 7: J367, Ui case: the region inside the annulus is where ^3(40102 — n^) > 0. 



Andrews' squeezing systi 




Andrews' squeezing syster 




-0.05 -0.04 -0.03 -0.02 -0.01 0.01 0.02 0.03 



Figure 8: Singular positions (according to J367, f/i) when bi = —0.02691, ai = 0.02, 02 = 
0.055 and 03, . . . , ay have the original values. The physical interpretation is as in Figure 

El 

6.1 Appendix 

The coefRcients fii The coefficients /i, . . . , /s in the context of T5 are 



/i = 4(05 - a^Yibi - 2biwi + &2 - + wi + w^) 

= 4(05 — 04)^16 — wp, 
/2 = 4(wi — 61) (04 — as) (—6^ + 2biWi — b\ + 262^2 — wl — + (i\ — + 2aia^ — a^) 

= 4(wi - 6i)(a4 - as) (03 - («4 - a^f - \b - wp). 



/3 


= bl 


- 2biWi + bl 


— 262^^2 + 262a4 — 262as + w'l + wl — 2w2a4^ + 2w2a5 


— a\ + a\ — 2a4a5 + ag 




= \b 


-wp + 2(62- 


- W2)(a4 - as) - ^3 + («4 - a^y, 




/4 


= b\ 


- 2biwi + bl 


— 262W2 — 262a4 + 262as + wl + wl + 2w2a4 — 2w2as 


— a^ + a\ — 2a4a5 + ag 




= \b 


-w\^-2{b2- 


- u^2)(a4 — as) — ag + (a4 — as)^. 




h 


= al 


— a\ + 2a4as 


~ ^5 ~ bl + 2biWi — &2 + 2&2'U^2 ~ w\ — wl 






= «3 


- (a4 - as)^ 







Expressed in angles, 



variables 


singularity 1 


singularity 2 


Cl 


-0.3621 


0.0127 


■'''1 


-0.9322 


0.9999 


C2 


0.1860 


0.1860 


S2 


0.9862 


-0.9826 


C3 


0.6364 


0.6364 


S3 


-0.7714 


-0.7714 


C4 








S4 


-1 


-1 


C5 


0.7714 


0.7714 


S5 


0.6364 


0.6364 


C6 








S6 


1 


1 


Cy 


-0.6364 


-0.6364 


S7 


0.7714 


0.7714 


se are 


Angles 


singularity 1 


singularity 2 


yi 


-1.9413 


1.5581 


y2 


1.3837 


-1.3837 


ys 


-0.8810 


-0.8810 


y.i 


1.5708 


1.5708 


Z/5 


0.6898 


0.6898 


1/6 


1.5708 


1.5708 


yi 


2.2606 


2.2606 



Table 2: The singularities of J367 type, original values apart from 61,01,02. The values 
are presented only with 4 decimals but were computed with 16 decimals. 

The coefficients dj, Iji The coefficients dj, Ij in the context of K2 ar e 

di = 2aia2(a| 2036153 - 20362C3 + h\^- h\) 
d2 — — 4ofo|s| 

ds — —of + 2a\a2 + a\a\ -\- 20^0361^3 — 20^036203 + a\h\ 
-\-a\hl — 02 + O20I + 20I0361S3 — 202O362C3 + albl + o|6| 



h = -20102(03 + 20361S3 - 20362C3 + bl + h\) 

I2 =20102(0353-^61) 

^3 = -(ascs - 62) (o? -al + al + 20361S3 - 20362C3 + bj + hi) 

I4 = 2O1O2S1S2 - (O3S3 + 6i)02C2 + (O3C3 - 62)0252 - (O3S3 + 6i)Oi. 




Figure 9: A singular configuration with rational Ci,Si,ai,bi. 



We can also simplify tlicsc expressions: 



do 


= |6|2 + 2a3(6iS3-&2C3) 


di 


= 2aia2(io 


d2 


= nin2 


ds 


— {a\ + a^do — {a\ — 


ni 


= (ai + 02)^ - do 


n2 


= (oi - 02)^ - do = 4aia2 - rii 


h 


— —di 


h 


= -{a3C3-b2){al- al + do) 


h 


= -{0,383 + hi){a2C2 + ai) + a2S2{a3C3 - 62 + 2aiSi) 


9i 


= —'ia\a\s\ + ni{Aaia2 — rii) 


92 


= 2aia2doC2 + (fli + a2)do — (a^ — alY 


93 


= -2aia2doSi + 20102(0353 + 61) - (0303 - 62) (oi - al + do) 




= {al — al)ci + U 



The coefficients r^: The coefficients Vj in the context of Li are 



ri = {aj + 02)1^1^ ~ 2biaja3S7 — 2biala3S7 + 2b2aja3C7 + 262020307 — {of — a^)^ + {aj + al)al 
r2 = 2ai{bia2 - 020357)^2 

r3 = 6^62 + b\a^c-j — 261620357 — 2bia\c-js-j + b\ + 3630307 + 620^ — 62O2 + 362O3C7 

+62O3S7 + ala^ci — a^asCi + a^cr 
r4 = (2aia2)siS2 + (-6102 + 020387)02 + (-62O2 - 020307)52 + (-6iai + 010357) 



References 

[ArpOl] T. Arponen. Regularization of constraint singularities in multibody systems. 
Multibody System Dynamics, 6(4):355-375, 2001. 

[BA94] E. Bayo and A. Avello. Singularity-free augmented Lagrangian algorithms for 
constrained multibody dynamics. Nonl. Dyn., 5:209-231, 1994. 

[CL092] D. Cox, J. Little, and D. O'Shea. Ideals, Varieties and Algorithms. Springer, 
Berlin, 1992. 

[EH95] E. Eich and M. Hanke. Regularization methods for constrained mechanical multi- 
body systems. Z. Angew. Math. Meek., 10:761-773, 1995. 

[Eis96] D. Eisenbud. Commutative Algebra, volume 150 of Graduate Texts in Mathemat- 
ics. Springer, 1996. corr. 2nd printing. 

[Gil78] D.R.A Giles. A comparison of three problem- oriented simulation programs for 
dynamic mechanical systems. PhD thesis. University of Waterloo, Waterloo, Ontario, 
1978. 

[GP02] G.-M. Greuel and G. Pfister. A Singular introduction to commutative algebra. 
Springer- Verlag, Berlin, 2002. With contributions by Olaf Bachmann, Christoph 
Lossen and Hans Schonemann, With 1 CD-ROM (Windows, Macintosh, and UNIX). 

[GPS05] G.-M. Greuel, G. Pfister, and H. Schonemann. Singular 3.0. A Computer Alge- 
bra System for Polynomial Computations, Centre for Computer Algebra, University 
of Kaiserslautern, 2005. http://www.singular.uni-kl.de 

[HW91] E. Hairer and G. Wanner. Solving Ordinary Differential Equations II: stiff and 
differential-algebraic problems, volume 14 of Computational Mathematics. Springer, 
1991. 

[LR05] B. Leimkuhler and S. Reich. Simulating Hamiltonian Dynamics. Cambridge Uni- 
versity Press, 2005. 

[ManSl] D.W. Manning. A computer technique for simulating dynamic multibody systems 
based on dynamic formalism. PhD thesis. University of Waterloo, Waterloo, Ontario, 
1981. 



[McCOO] J.M. McCarthy. Geometric Design of Linkages, volume 11 of Interdisciplinary 
Applied Mathematics. Springer, New York, 2000. 



[MI03] F. Mazzia and F. lavernaro. Test set for initial value problem solvers. 
Release 2.2, Department of Mathematics, University of Bari, August 2003. 
[http: / / pitagora.dm.uniba.it / ^testset / , 

[Nor 76] D. Northcott. Finite Free Resolutions, volume 71 of Cambridge Tracts in Math- 
ematics. Cambridge Univ. Press, 1976. 

[RS88] R. Roberson and R. Schwertassek. Dynamics of Multibody Systems. Springer, 
1988. 

[Sch90] W. Schiehlen, editor. Multibody Systems Handbook. Springer, Berlin, 1990. 



